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ABSTRACT 

We use a semi-analytic model of Lya em.itters (LAEs) to constrain the reionization 
history. By considering two physically motivated scenarios in which reionization ends 
either early (ERM, Zi « 7) or late (LRM, Zi « 6), we fix the global value of the IGM 
neutral fraction (e.g. xhi — 3x 10"'', 0.15 at z = 6.56 for the ERM and LRM, respec- 
tively) leaving only the star formation efficiency and the effective escape fraction of 
Lya photons as free parameters. The ERM fits the observed LAE luminosity function 
(LF) at z = 5.7 and 6.56 requiring no redshift evolution or mass dependence of the star 
formation efficiency, and LAE star formation rates (SFR) of 3 < M*/M0yr~^ < 103, 
contributing « 8% of the cosmic SFR density at z = 5.7. The LRM requires a phys- 
ically uncomfortable drop of « 4.5 times in the SFR of the emitters from z — 6.5 
to 5.7. Thus, the data seem to imply that the Universe was already highly ionized 
ai z = 6.56. The mass-dependent Lya transmissivity is 0.36 ^ Ta ^ 0.51 (ERM) and 
Ta ^ 0.26 (LRM) at z = 6.56. The LF data at z = 4.5 imply an extra Lya line damp- 
ing factor of w 0.25 possibly due to dust; the presence of a (clumpy) dust component 
with E{B — V)^ 0.28 is also required to reproduce the observed large Lya equivalent 
widths at the same redshift. Additional useful information can be extracted from the 
line profile (weighted) skewness, found to be Sw = 10 — 17 A for the two reioniza- 
tion models, which shows an interesting L^ — Xhi anti-correlation, holding under the 
model assumptions. The shortcomings of the model and strategies to overcome them 
are discussed. 

Key words: linc:profiles - galaxies:high redshift - luminosity function - intergalactic 
medium - cosmology:theory 



1 INTRODUCTION al. 2007, Spergel et al. 2007) and high redshift galaxy surveys 

(Bouwens et al. 2006, Stark et al. 2007). 



The Epoch of Reionization (EoR) marks the second ma- 
jor change in the ionization state of the universe after re- 
combination and is directly linked to structure formation. 
Reionization begins when the first structures form within 
dark matter halos and emit neutral hydrogen ionizing pho- 
tons. In addition to changing the ionization state, these first 
structures also affect subsequent structure formation due to 
various radiative, mechanical and chemical feedback effects. 
Thus, to probe reionization, one needs an excellent under- 
standing of initial density perturbations and their growth, as 
well as simulations that can trace the evolution of structure 



The emerging picture (Choudhury & Ferrara 2007) is 
one in which hydrogen reionization is an extended process 
starting at z « 15 and being 90% complete hy z = 8. Reion- 
ization is initially driven by metal-free stars in low mass 
(A/ < 10* M0) halos; the conditions for the formation of 
these objects are soon erased by the combined action of 
chemical and radiative feedbacks ai z < 10. 

Given the many assumptions necessarily made by reion- 
ization models, the above scenario needs constant confronta- 
tion with freshly acquired data sets. In this sense, it has been 



r . suggested (Malhotra & Rhoads 2004, 2005; Santos 2004; 

„ „ . , . , „ „ . . , . 11-, Haiman & Cen 2005; Mesinger, Haiman & Cen 2004; Di- 

One ot the major challenges ot reionization models IS to ,,, , n tt t-n 1 

, ,,,.,, , , r ,1 -111 1 ikstra, Wyithe & Haiman 2007; Diikstra, Lidz & Wyithe 

be able to simultaneously account tor the considerable, and 1„„„ \r ■ » r-, n ,,^^.,n 1 , r , • / -, 

r, ,, n- ,• . r 1 . i . , i 2007; Mesmger & lurlanetto 2007) that a class oi high red- 

orten apparently connictmg, amount ot data accumulated by , r i i -r » i i t^ 

, .,. r^rir^ i ,■ i- . /T-, , shift galaxies, the Lyman Alpha Emitters (LAEs) can be 

experiments exploiting QSO absorption line spectra (Ian et ^ ' 

1 nr,n^\ ■ ■ 11 1 1- i- /n i Suitably used to put additional constraints on the reiomza- 

al. 2006), cosmic microwave background radiation (Page et . , . , t ■, -, , , ,■■,-, 

tion history: the Lyman break and the strength, width and 

asymmetry of the observed Lya line make the detection of 

* E-mail: dayal@sissa.it (PD) LAEs unambiguous. The strength of the method is based 
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on the sensitivity of Lya photons to even tiny amounts of 
H I in the intergalactic medium (IGM). At redshifts 2 ~ 5, 
the optical depth to Lya photons is very large. Let e be the 
electron charge, / the oscillator strength (0.4162), A^ the 
wavelength of Lya in its rest frame (1216 A), me the elec- 
tron mass, c the speed of light, H[z) the Hubble parameter 
at the required redshift, uhi the global neutral hydrogen 
density and uh the global mean hydrogen density at that 
redshift. Note that uhi ~ XHinn where xhi is the fraction 
of neutral hydrogen at the redshift under consideration. Fur- 
ther, Qb represents the baryonic density parameter and Qrn 
is the total (baryonic -I- dark) matter density parameter of 
the universe, Sim = Sib + ^dm- Then, 

ne fXa riHi 
meCH(z) Uh 



where 



7re f\a 
m^cH(z 



■UH = 1.76 X 10^/i"^S7„^/^ 



Qth i 1 + z 



0.022 



8 



3/2 



Hence, even a H i fraction of 10~^ can lead to a significant 
attenuation of the Lya line. The observed (i.e. transmitted) 
Lya luminosity. La , can then be used to infer the ionization 
state of the IGM at redshifts close to those of the emit- 
ter and hence to reconstruct, at least piecewise, the cosmic 
reionization history. 

This simple picture is complicated by a number of im- 
portant physical effects. First of all, Lya photons from the 
stars have to propagate through and escape from the inter- 
stellar medium of the LAE. During their travel they are mul- 
tiply scattered by H i atoms (thus being either removed from 
or added to the line of sight [LOS]) and possibly absorbed 
by dust grains (Neufeld 1991; Tasitsiomi 2005; Hansen & Oh 
2006; Finkelstein et al. 2007). These processes modify both 
the emerging Lya luminosity and the shape and equivalent 
width of the line. Second, the ionizing radiation from the 
same stars builds regions of ionized IGM around the emit- 
ters, whose size depends on the star formation rate, age, es- 
cape of ionizing photons from the galaxy and the stellar Ini- 
tial Mass Function (IMF; the case of very massive stars has 
been explored, for example, by Dijkstra & Wyithe 2007). As 
a result, the flux redwards of the Lya line can escape, attenu- 
ated only by the red damping wing of the Gunn-Peterson ab- 
sorption (Miralda-Escude 1998; Madau & Rees 2000). To a 
first approximation, the spatial scale imposed by the Gunn- 
Peterson damping wing on the size of the H 11 region corre- 
sponds to a redshift separation of Az ~ 0.01, i.e. about 200 
kpc (physical) at z = 10. The effects of the damping wing 
fade away if the emitter is powerful enough to create a large 
enough H 11 region and/or if the universe is already reion- 
ized when the emitter turns on. Alternatively, one would 
observe the damping wing if there were even a small frac- 
tion of neutral hydrogen left inside the sphere and/or if a 
H I cloud is present along the LOS to the source. 

All the above effects combine to shape the observed 
LAE Luminosity Function (LF), which has been now mea- 
sured (Rhoads et al. 2000; Taniguchi et al. 2005; Shimasaku 
et al. 2006; lye et al. 2006; Kashikawa et al. 2006; Murayama 
et al. 2007; Ota et al. 2007; Dawson et al. 2007) with different 
degrees of accuracy up to « « 7. Such tremendous progress 
has been made possible by the increase of survey fields and 
available samples. The current observational situation can 



be summarized as follows. All studies seem to converge to- 
ward the conclusion that there is very little indication of 
evolution of the LF moving from z — 3 to z — 5.7. Beyond 
that epoch there seems to be evidence of a decline in the LF, 
with L^, a,t z = 6.6 being about 50% of that a,t z — 5.7. Such 
a high luminosity steepening of the LF can be produced by 
a number of different physical effects. A rapid evolution of 
the IGM ionization state can be invoked (Kashikawa et al. 
2006) if the overlapping phase of reionization ended around 
z — 6; however, the net effect of reionization on the observed 
Lya luminosity of the most luminous (and presumably mas- 
sive) LAEs is unclear. If these objects are expected to live 
in more dense and hence more neutral environments, they 
are also more heavily clustered (McQuinn et al. 2007). The 
two effects might not change appreciably the size of their 
H II regions. Alternatively, the observed evolution could be 
simply a result of the evolution of the mass function of dark 
matter halos housing the LAEs (Dijkstra, Wyithe & Haiman 
2007; Dijkstra, Lidz & Wyithe 2007). Finally, extinction due 
to dust, which is expected to be more prominent in actively 
star forming galaxies, may act as a sink for Lya photons in 
the most luminous LAEs. 

As of now, it is difficult to firmly assess which of these 
explanations is more robust. Fortunately, other aspects of 
the data, such as the line shape and equivalent width, might 
allow one to make progresses. Here we try to assess to what 
extent the reionization history can affect the shape of the 
LF and the observed properties of individual LAEs. Our ap- 
proach is similar in spirit to some of those mentioned above, 
but it has the strength of being based on reionization models 
that simultaneously account for all the available data beyond 
LAEs including Lya/Ly/3 Gunn-Peterson opacity, electron 
scattering optical depth, Lyman Limit Systems, cosmic SFR 
history and the number density of high-redshift sourcesQ 



2 THE MODEL 

In this section we describe the physical features of the model 
we have developed to derive the various properties of LAEs 
which will then be compared with observations. Several 
steps are required in order to carry out this task which are 
described in detail in the following. These include the use 
of the Sheth-Tormen mass function to obtain the redshift 
dependence of the number density of dark matter halos, the 
star formation prescriptions required to build the luminosity 
function, the production rate of H i ionizing photons and the 
intrinsic Lya luminosity, the size of the Stromgren sphere 
built by LAEs and the H i density profile within it and in 
the general IGM, for which we use a previously developed 
reionization model. 



^ Throughout the paper, we use the best-fit cosmological param- 
eters from the 3-year WMAP data (Spergel et al. 2007), i.e., a 
fiat universe with {Um, ^A, Hj/i^, /i)=(0.24, 0.76, 0.022, 0.72). 
The parameters defining the linear dark matter power spectrum 
are erg = 0.82, Us = 0.95, dns/dlnk = 0. We use a value of rrg 
much higher that quoted from WMAP3 (0.76) as the combination 
of WMAP3 and SDSS data give ag ~ 0.78 (0.86) for low (high) 
resolution Lyo forest data (Viel et al. 2006). Mpc is comoving 
unless otherwise specified. 
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2.1 The mass function 

We start with the well known Sheth-Tormen mass function, 
ISheth fc Tormenl (|l999l ) , which is used to calculate the num- 
ber density of dark matter halos of mass between M and 
M + dM at any redshift z, represented by n{M, z)dM, as 



niM,z)dM = A[l + ^],l-^—. 



1 



2 p dry' 



-i/'^/2 



rfM, 



(1) 



where v' — ^/av. 

In eq[T] A, a and q are modifications to th e original 
Press-Schechter mass function. IPress fc Schechteij (11974), to 
make it agree better with simulations. Here, A « 0.322, 
q = 0.3 and a = 0.707. 

As in the Press-Schechter mass function, 



/* 1 ^6 ,^ 



(4) 



Using the popu lation synthesis code Starburst99 
(JLeitherer et al.ll 19991 ) we obtain the hydrogen ionizing pho- 
ton rate, Q, emitted by galaxies having a given SFR, assum- 
ing a metallicity Z = O.OSZ© . Determining the metallicity of 
the LAEs proves very challenging, as for most of the cases, 
only the Lya line can be detected from these objects. To 
guess their metallicity, we use the results from studies of 
LBGs (Lyman Break Galaxies) and DLA (Damped Lya) 
systems, which indicate values of 0.05 — G.IQZq, which jus- 
tifies our assumption, iPettini (2003). We use a Kroupa IMF 
with a slope of 1.3 between 0.1 and O.SM© and 2.35 between 
0.5 and lOOM©. Using the fact that Q scales linearly with 
SFR, we can calculate Q for the desired SFR. 



D{z) 



D{z)g{M)' 
ffW/[5(0)(l- 



nA + (l+!^m/2)(l+nA/70)]" 



Here, 5c{= 1.69) is the critical overdensity for spherical col- 
lapse and D{z) is the growth factor for linear fluctuations, 
ICarroll. Press fc Turneij (|l992l ). Further, the variance of the 
mass M contained in a radius R is given by 



-^ (-R) = ^ / k'P[k)W\kR) — . 



(2) 



In eqU W{kR) = 3(sin(fc_R) - kRcos{kR)) is the window 
function that represents the Fourier transform of a spherical 
top hat fllter of radius R, P{k) = Apk"T^(k) is the power 
spectrum of the density fluctuations, extrapolated to z = 
using linear theory where Ap is the amplitude of the density 
fluctuations calculated by normalizing a{M) to ag which 
represents the variance of mass in a sphere of size 8h~^ Mpc 
at z = 0. The term T{k) is a transfer funct ion which rep- 
resen ts differential growth from early times (jBardeen et al.l 
1 19861 ). 



T(fc) = 



0A3q-^ ln(l + 2.34g) 



[1 + 3.89g + (16.1g)2 + (5.46g)3 + (6.71g)4]i/4 ' 



(3) 



where q = k{flrnh'^)~^ ■ 

Once the mass function is obtained, a SFR recipe (Sec 
2.2) is used to obtain the intrinsic Lya luminosity for any 
halo on the mass function in Sec 2.3, thereby providing the 
intrinsic Lya luminosity function. The attenuation of the 
intrinsic Lya luminosity by the IGM, as calculated in Sec 
2.4, then allows the mass function to be translated into the 
observed Lya luminosity function. 



2.2 The ionizing photon rate 

The baryonic mass, Mt, contained within a halo of mass Mh 
can be expressed as 



Mb 



^6 



-Mh. 



We assume that a fraction /, of this baryonic matter forms 
stars over a timescale i, = edc^H , where edc is the duty cycle 
and tn is the Hubble time at 2: = 0. Thus, we can write the 
star formation rate (SFR) as 



2.3 Intrinsic Lya line 

Star formation in LAEs produces photons with energy > 
1 Ryd. These photons ionize the interstellar H i , leading 
to the formation of free electrons and protons inside the 
emitter. Due to the high density of the ISM, these then 
recombine on the recombination time scale, giving rise to a 
Lya emission line. 

Let fesc be the fraction of H i ionizing photons that 
escape the galaxy without causing any ionizations, fa the 
fraction of Lya photons that escape the galaxy without be- 
ing destroyed by dust, Va be the frequency of Lya in the rest 
frame of the galaxy (1216 A) and h be the Planck constant. 
Then, the intrinsic Lya luminosity, I/™*, from the galaxy 
can be expressed as 



tQ(1 - fesc)fahVa. 



(5) 



It has been calculated that there is a two-thirds probability 
of the recombination leading to a Lya line and a one-third 
probability of obtaining photons of frequencies different from 
the Lya (Ostorbrock 1989). This gives rise to the factor of 
two-thirds in eq[S] For (1 — /eac)/a = 1, the intrinsic Lya 
luminosity and the SFR are related by the following 

rint r, on , 1 n42 -1 SFR 

La = 2.80 X 10 erg s — 

Moyr-i 

Modeling the Lya line to be Doppler broadened, the 
complete line profile is 



2 1 

S y'-r^Ai'd 



exp 



-{,'-^a.r/A,.j 



(6) 



where Aj^d = {vc/c)h'a, Vc is the rotation velocity of the 
galaxy and c is the speed of light. 

The minimum rotation velocity of the galaxy would be 
equal to the rotation velocity of the host halo, Vh- How- 
ever, for more quiescent star formation, for realistic halo 
and disk properties, Vc can have values bet ween Vh and 
2vh (|Mo. Mao fc Whitdl 19981 : ICole et al.ll2000l ). We use the 
middle value between these limits in our model, so that 
Vc = l-5vh- To illustrate, as Mh increases from 10^'' to 
10^'^ Mq, Vc increases from 102 to 475 km s~^ at z ~ 6.6. 

We calculate the velocity of the halo assuming that the 
collapsed region has an overdensity of roughly 200 times the 
mean cosmic density contained in a radius r2oo- Then, v^, 
the velocity at r2oo is expressed as 
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vl{z) 



GMh 

r200 



= GMh 



GMh 



1/3 



(7) 



where uhi, rir, 



are the number density of neutral hy- 



where Q,m and H are the density and Hubble parameters, 
respectively, at the redshift of the emitter. 

To summarize, the intrinsic Lya luminosity depends 
upon: the ionization rate Q, the escape fraction of 
H I ionizing photons fesc, the escape fraction of Lya pho- 
tons fa and the rotation velocity of the galaxy Vc . In turn, 
Q depends on the SFR (which is a function of halo mass), 
the metallicity Z, and the age of the emitter t« , chosen such 
that the number of ionizing photons emitted per second set- 
tles to a constant value. 



2.4 Observed Lya line 

The intrinsic Lya line is attenuated by the neutral hydrogen 
present in the IGM along the line of sight toward the emitter. 
In this section we compute the neutral hydrogen distribution 
and the attenuation caused by it. 



2.4-1 Global xh I calculation 

We use the global value of the H i fraction 
Y-ffj = UHilnH resulting fr om the modeling by 



iGallerani. Choudhury fc Ferraral (|2006l ). further refined 
in iGallerani et al.l ((20071) • The main features of the model 
are summarized here. Mildly non-linear density fluctuations 
giving rise to spectral absorption features in the Inter- 
galactic medium (IGM) are described by a Log-Normal 
distribution. This has been shown to fit the observed 
probability distribution function of the transmitted flux 
between redshifts 1.7 and 5.8 bv iBecker. Ranch fc Sargent) 
((20071). For a given IGM equation of state, this being the 
temperature-density relation, the mean global H i fraction 
(Xhi) can be computed from photoionization equilibrium 
as a function of baryonic over-density (A = p/p) and pho- 
toionization rate (Fb) due to the ultra-violet background 
radiation field. These quantities must be d etermined from 
a com bination of theory and observations. iGallerani et al.l 
((20071) included two types of ultraviolet photons: from 
QSOs and Pop II stars. The free parameters in their model 
were (i) the SFR efficiency (/*) and (ii) the escape fraction 
of ionizing photons from the galaxy (fesc)- These were 
calibrated to match the redshift evolution of Lyman-limit 
systems, Lya and Ly/3 optical depths, electron scattering 
optical depth, cosmic SFR history and number density 
of high redshift sources. The following reionization sce- 
narios provide a good fit to observational data: (i) Early 
Reionization Model (ERM), in which reionization ends at 
Zi = 7, (/* = 0.1, fesc = 0.07), (ii) Late Reionization Model 
(LRM), where reionization ends at Zi = 6, (/* — 0.08, 
fesc = 0.04). 



2.4- S Neutral hydrogen profile 

The IGM is approximately in local photoionization equi- 
librium. Under such conditions ionizations are balanced by 
recombinations, 



drogen, protons and electrons respectively, as is the hydro- 
gen Case B recombination coefficient and F_b is the ioniza- 
tion rate due to the background. As mentioned in Sec l2.4.ll 
in this work we take advantage of the results presented by 
iGallerani et al.l ((2003). Once that xhi is fixed to the their 
valueo the photoionization rate contributed by the ionizing 
background light produced by quasars and galaxies is given 
by: 



rB = 



XHI 



-uholb- 



(9) 



Moreover, the radiation from stars inside the galaxy ionizes 
the region surrounding the emitter, the so-called Stromgren 
sphere. 

The evolution of the Str omgren sphere is given 
by the following relation , ((Shapiro fc Girouxl il987i : 
iMadau. Haardt fcReeJi 19991 ) 



-J- - 3H{z)Vi = 

dt Uhi 



Vi 



(10) 



where, Vi is the proper volume of the Stromgren sphere, and 
tree = [1.1 7aBnp]~ is the volu me averaged recombination 
timescale (JMadau fc Reesli200(]| ). The proper radius Ri = 
(3V7/47r)^'^, identifies a redshift interval Az between the 
emitter and the edge of the Stromgren sphere, given by the 
following: 



Az = 100{nmh'^y^^{l + zf^Rilc. 



(11) 



Though this equation is not strictly valid at z ~ 0, it is a 
good approximation at the high redshifts we are interested 
in (z J5 4.5). If Ze is the redshift of the emitter, for redshifts 
lower than the Stromgren s phere redshift, i.e . . Zs = Ze — Az, 
we use the xhi value from iGallerani et al.i (20071). Within 
the Stromgren sphere, to Fs we add the LAE photoioniza- 
tion rate F^: 



-sir 



Aur 



ro-i 



he 



dA, 



(12) 



where Lx is the specific ionizing luminosity of the emitter 
(in erg s~^A~^), \l is the Lyman limit wavelength (912 A) 
and gl is the hydrogen photoionization cross-section. Thus, 
inside the ionized region, xhi is computed as following: 



XHi(r) = 



2nHCtB + r(r) ± ^r2(r)+4nj:/asF(r) 
2nHOiB 



(13) 



where T{r) = F_B(r) -fFs. The solution must be chosen such 
that Xhi < 1, which only happens for a negative sign before 
the square root. At the edge of the Stromgren sphere, we 
force XHi{r) to attain the global value in the IGM. 



2.4-3 Lya optical depth and transmitted flux 

The transmitted Lya luminosity is Ta ~ e~'^° where t„ is 
the optical depth to Lya photons. Assuming that reioniza- 
tion completes Skt z = Zi, Tc, can be calculated as 

P' dr 

Ta{Vobs) = / cr{Uobs)nHi{z)--dz, 
/ dz 



nniT'B ~ UeUpaB, 



(8) 



We assume a homogenous and isotropic IGM density field. 
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ao(f>{uobs)nHi{z) — dz, 



where a is the total absorption cross-section, ctq = 
7re^//(mec) and 4> is the Voigt profile. 

For regions of low H i density, the natural line broad- 
ening is not very important and the Voigt profile can be 
approximated by the Gaussian core: 



V^Aud 



exp 



-{•^ob.-''o,.^r/^''i 



(14) 



In ea ll4l Va^,- — c/[Xa{l + Zr)] is the local Lya frequency at 
a distance r from the emitter and Vots = c/Xobs- Further, 
Aud = b/Xa, where 6 — W2kT/rnH is the Doppler width 
parameter, niH is the hydrogen mass, k is the Bol tzmann 
constant and T = W^K is the IGM temperatur e (jSantosI 
I2OO4I : ISchave et al]|2000l : iBolton fc Haehneltll2007l ) . 

For regions of high H i density, we take into account 
the Lorentzian damping wing of the Voigt profile. Thus, for 



wavelengths outside the Gaussian cor e, i.e. for | 1/ - 



Aud, we assume the following profile (|Peebleslll993l ): 



^Lorentz — 



'i-K'^{vobs — y 
where A = 8.25 x 10** s" 



\6' 



1^ 



(15) 



,,,)' + {l^^l^)[v,bslv^ 

^ is the decay constant for the Lya 



resonance. 



3 BASIC DEPENDENCIES 

By using the model described in the previous section, we 
can compute the observed Lya line profile; 



^ — ^a r ^nt rp r %'. 



(16) 



The Lya optical depth depends on three quantities: the 
star formation rate (which fixes the value of Q), the ion- 
ized region radius, and the global neutral fraction: r^ = 
rc,{M-t , Ri , Xh i) ■ Once these three parameters are given, 
the transmissivity is uniquely determined. Notice that Ri — 
Riifesctt,, Mi,,XHi)- If instead we are interested in the ob- 
served Lya luminosity, a fourth parameter needs to be spec- 
ified, the "effective" Lya photon escape fraction 



fe. 



Jescjja, 



(17) 



which expresses the physical fact that the condition to ob- 
served Lya photons is that some ionizing photons are ab- 
sorbed within the galaxy and only a fraction fa of produced 
Lya photons can escape to infinity. Note that fesc,a does 
not affect the transmissivity as both the intrinsic and the 
observed luminosity depend on it and therefore it factors 
out. A full exploration of the physical effects of the param- 
eters on the observed luminosity. La, can be performed by 
varying only the parameters M^,,RI,XHI and fssc,a- The 
effects of other parameters (as, for example, metallicity, Z) 
can be estimated by simple scaling of the results below. 

To understand the impact of each of the three relevant 
quantities on La we have selected a fiducial case with pa- 
rameters broadly similar to those we inferred under realistic 
(i.e. observationally derived) conditions for LAEs and allow 
them to vary in isolation taking three different values. We 
therefore considered 1 x fiducial -1-4x3 = 13 different cases 
shown in Fig.l and summarized in detail in Table 1. 



Table 1. Parameters of the fiducial model as well as for the 
diff'erent cases plotted in Fig[T] For all cases, the halo mass is 
10^^'^ Mq. Dashes indicate that fiducial model values have been 
used. 



Model 



Ah 
[Mgyr- 



XHI 



Ri 

[pMpc] 



Ta 



Fiducial 


27 


si 


81 


s2 


54 


s3 


13.5 



0.35 



0.01 



5.95 



0.47 
0.52 
0.50 
0.44 



fl 


- 


0.9 


— 


— 


0.47 


12 


- 


0.1 


— 


— 


0.47 


13 


- 


0.03 


- 


- 


0.47 


rl 


- 


- 


— 


2.97 


0.44 


r2 


- 


- 


- 


1.48 


0.37 


r3 


- 


- 


- 


0.74 


0.27 


cl 


_ 


_ 


3 X 10-4 


— 


0.49 


c2 


- 


- 


0.05 


— 


0.42 


c3 


- 


- 


0.15 


- 


0.32 



3.1 Star formation rate 

The ionizing photon rate, Q, of the emitter is directly pro- 
portional to its SFR. As a result, a larger SFR results in 
(a) an increase of I/^', (b) a larger ionized region around 
the LAE, (c) a lower value of xhi at each point within the 
Stromgren sphere (see eqs. 12-13). The net effect is that as 
SFR increases, the transmission of a stronger Lya line in- 
creases due to decreased damping by both the Gaussian core 
and the red damping wing. This is shown in panel (a) of 
FigUl For the fiducial case we find that 47% of the intrin- 
sic Lya luminosity is transmitted; this value increases with 
SFR, reaching 52% when M* = SIMqjt'^ , as seen from 
TabU] 



3.2 Effective Lya photon escape fraction 

The effective Lya photon escape fraction fe3c,a scales both 
Lq"* and La equally, without changing either the size of the 
Stromgren sphere or the H i profile within it. The fraction 
of Lya luminosity transmitted is hence, the same in all the 
cases. The variation of La with fssc,a is shown in panel (b) 
of FiglU 



3.3 Ionized region radius 

As the ionized region becomes larger, due to a more robust 
input on ionizing photons from the source, the Lya photons 
reach the edge of the sphere more redshifted. Hence, the 
H I outside the ionized bubble is less effective in attenuating 
the flux. The size of the ionized region radius is therefore 
very important for LAEs in regions of high H i density and 
loses importance as the H i density decreases. We show the 
variation of La with ijjj in panel (c) of Fig[T] from which we 
can readily appreciate that as Rj increases (at a fixed xhi 



Rj is in physical Mpc (pMpc) 
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Figure 1. Effect of varying (a) SFR (b) fesc,a (c) Ri and (d) XHI on La. Refer to Tab[T]for the parameters used for eacli of tiie lines 
in this plot. The dashed vertical line shows the wavelength of the redshifted (emission redshift z = 6.56) Ly« line. 



and SFR), a larger fraction of the line is transmitted due to 
the aforementioned effect. As, to a good approximation. 



Ri ex 



^ Jescf^* 

XHinH 



1/3 



(18) 



for a fixed value of Q (SFR) and xhi, Ri can vary either 
due to t, or /esc- These two parameters play a qualitatively 
different role. While the age variation can be embedded in 
a variation of Ri only, changing the value of the escape 
fraction also affects L^"* (see eq[6)l giving rise to a physically 
interesting effect. In FigO for illustration purposes, we fix 
AL = 27M0yr~\ U = 10® yr, fa = l and study the effect 
of fesc on La for different values of xhi- 

The observed Lya luminosity decreases monotonically 
with /esc for low values of xhi (< 0.01), just mirroring the 
decreasing value of the intrinsic Lya line. Here, the fact 
that the size of the Stromgren sphere built increases with 
increasing fesc has no effect on La simply because the H i 
density is too low to cause (red) damping wing absorption, 
irrespective of the size of the ionized region. For xhi ^ 0.01, 



the La trend with f^sc in not monotonic anymore (see also 
Santos, 2004). For example, for xm ~ 0-15 the observed 
Lya luminosity reaches a maximum at f^sc ~ 0.5. This can 
be explained by the following: for low (< 0.5) f^sc values, 
as /esc increases, the ionized volume increases, thus leading 
to larger transmission. When La reaches its maximum (for 
fesc ~ 0.5, in our example), a further fcsc increase reduces 
the observed Lya luminosity, as a consequence of the de- 
creasing value of L^"*. This highlights the fact that while 
for low values of Xm^ fesc affects the observed Lya only 
through the intrinsic Lya line, for high values of xhIj the 
effect of fesc on the Stromgren sphere size becomes consid- 
erably important. 

3.4 Neutral hydrogen fraction 

In panel (d) of Fig[T] we study the effect of different xhi 
values on the Lya line. It can be seen from Tab. [l] that 
the Lya line is quite damped {Ta ~ 0.32) for high values 
of Xhi {~ 0.15). As the value of xhi decreases, the effect 
of both the Gaussian core and the red damping wing start 
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candidates; 73 of which were confirmed using DEIMOS on 
KECK II and the Low Resolution Imaging Spectrograph 
(LRIS). 

IShimasaku et alj (|2006l ) identified 89 LAE candidates 
in the Subaru Deep Field (SDF) at z ~ 5.7 by using the 
8.2m Subaru Telescope and the following selection criteria: 
(a) i' - NB816 ^ 1; (b) NB816 ^ 26. By using the Faint 
Object Camera and Spectrograph (FOCAS) on Subaru and 
DEI MOS, 28 candidates w ere confirmed as LAEs. 

iTaniguchi et al.l l|2005r ) detected 58 possible LAEs using 
Subaru at z ~ 6.5 and obtained the spectra for 20 of them 
using the FOCAS. They found that only 9 of the above 
objects showed sharp cut-off at the Lya wavelength, narrow 
line widths and asymmetric profiles, thus being confirmed as 
L AEs at z ^ 6.5. These included the two LAEs discovered 
bv lKodaira et al.1 l|2005h at z = 6.541 and 6.578. Using the 
same selection criterion and instruments as .Taniguchi et al.l 
(|2005l ) and including the LAEs confirmed using the Keck 
II DEIMOS spectrograph, iKashikawa et all (|2006h added 8 
more LAEs at 6.5 to this list. Thus, the Subaru observations 
have a total of 17 confirmed LAEs at 2: ~ 6.5. 



Figure 2. Dependence of La on fesc for different values of XHI- 
Adopted parameters are M* = 27M0yr~l, U = 10® yr, fa = 1. 
The solid line shows the intrinsic Lyo luminosity. Curves with 
symbols refer to different values of Xm = 0.15,0.01, 10~^, 3 X 
10"'* from bottom to top, respectively. 



reducing, allowing more of the line to be transmitted. For 
Xhi = 3 X 10"*, most of the line redwards of the Lya wave- 
length escapes without being damped. This occurs because 
the emitter is able to (a) strongly ionize the H i within the 
Stromgren sphere (already ionized to a large extent even out- 
side it) even further, and (b) build a large Stromgren sphere 
such that the Lya line is not affected by the damping wing 
of the H I outside. 

We remind the reader that La = 

La{Mi,, fcac,a,Ri,XHi)- For a continuous star forma- 
tion mode, the luminosity of the source becomes rapidly 
independent of age (typically after 100 Myr) ; if, in addition, 
we ad opt the values of xhi obtained from ICallerani et al.l 
l|2007l) by matching the experimental data, we are left 
with two free parameters, M* and fesc,a- Recalling that 
Mi, oc ft/tdc, the free parameters in our model reduce to 
(a) /*/edc and (b) fesc,a- 



4 COMPARISON WITH OBSERVATIONS 

In this section we compare the results obtained from our 
model to observations of the LAE LF, the UV LF, the 
line profile asymmetries, the equivalent widths and the cos- 
mic SFR density. In particular, we would like to assess to 
what extent the study of these quantities for LAEs can 
be used to discriminate between the early (ERM) and late 
(LRM) reioniza t ion sc enarios, as deduced from the study of 
ICallerani et al.l |20o3), summarized in Sec l2.4.1l 



4.1 Available data 



iDawson et al.l (|2007|) conducted the Large Area Lya 
(LALA) survey to look for LAEs at 2: ~ 4.5 and found 97 



4.2 Lya Luminosity function 

As a first remark, it is useful to point out that if the LF evo- 
lution were to result purely from the evolution of the dark 
matter halos predicted by hierarchical structure formation, 
one would expect the comoving number density of luminous 
objects to increase with decreasing redshift. Although data 
errors are still large, it must be noted that instead there is 
an indication that there is no evolution of the Lya LF be- 
tween z ~ 3 - 6 (Dawson et al., 2007; Ouchi et al., 2007). 
Obviously, a number of different effects could produce this 
non-monotonic trend, a few examples being, SFR evolution, 
redshift dependent escape fractions and dust extinction, as 
we discuss in the following. In Fig|31 we plot the cumula- 
tive LFs at z = 4.5, 5.7 and 6.56 together with our best fit 
results. We now discuss the predictions of ERM and LRM 
separately. 

The ERM predicts an evolution of the hydrogen neutral 
fraction such that xhi = 1.3 x 10^^,8.6 x 10"^, 3 x 10"" 
for z — 4.5, 5.7 and 6.56 respectively. Interestingly, a very 
good fit to the data can be obtained for the two highest 
redshifts with a single value of the star formation efficiency 
parameter /./tdc ~ 3.5, thus implying that the SFR for any 
given halo mass is not very much dependent on redshift. 
While a reasonable fit to the data at z — 5.7 and z — 6.56 
is obtained for a single value of fesc,a ~ 0.3; a better fit is 
obtained by allowing for a 40% increase of fesc,a towards 
larger masses. The typical LAE dark matter halo masses 
corresponding to the observed luminosities are in the range 
Mh = 10^°-^""°Mo at z = 6.56; at the same redshift the 
star formation ranges from 2 to 43 M0yr~^. 

The data at z = 4.5 instead pose a challenge to the 
model because, assuming non-evolving values of ft/edc = 3.5 
and fesc,a, the observed number density of luminous objects 
is lower than that predicted by the evolution of the theoret- 
ical LF. Given the relative constancy of the star formation 
efficiency and of the effective Lya photons escape fraction 
noted for the two highest redshifts considered, the most nat- 
ural explanation is in terms of increasing dust extinction. 
To reconcile the prediction with the data at z — 4.5 we 



p. Dayal, A. Ferrara & S. Gallerani 

-2.5 



-3 = 



-3.5 






A 



-4.5 



O 



-5.5 




42 



42.5 



43 



43.5 



Log L„ [erg s-'] 



Figure 3. Cumulative LAE Luminosity Function for the early reionization model (ERM). Points represent the data at three different 
redshifts: z = 4.5 Dawson et al. (2007) (squares), z = 5.7 Shimasaku et al. (2006) (circles), z = 6.56 Kashikawa et al. (2006) with 
downward (upward) triangles showing the upper (lower) limits. Lines refer to model predictions at the same redshifts: z = 4.5 (dashed), 
z = 5.7 (dot-dashed), z = 6.56 (solid). 



then require that the Lya line suffers an additional damp- 
ing due to the presence of dust; which we find to be equal to 
1/4.0 = 0.25, i.e fa (and hence fesc,a) decreases by a factor 
of 4. A strong increase of the dust content inside galaxies 
is expected on cosmic time scales larger than 1 Gyr (cor- 
responding to z ;S 5) when evolved stars rather than core- 
collapse supernovae become the primary dust factories. Such 
a hypothesis needs to be checked carefully, as the dust would 
not only affect the Lya line but also the continuum emission, 
finally affecting the equivalent width of the line. We will dis- 
cuss these effects of dust in Sec. 4. 3 and 14.51 Hence, it seems 
that overall, a model in which reionization was completed 
relatively early {zi — 7) matches the data quite well. 

The LRM has a much slower reionization history, as is 
clear from the values of xhi = 1-4 x 10~^, 1.3 x 10"*, 0.15 
for z — 4.5, 5.7 and 6.56 respectively. At the lowest redshifts 
{z — 4.5 and 5.7) this model requires exactly the same value 
f*/^dc = 3.5 as the ERM. This does not come as a sur- 
prise of course, as xhi is so small at these epochs in both 



the ERM and the LRM that the observed Lya luminosity 
is unaffected. However, as xhi is much larger ai z — 6.56 
in the LRM as compared to the ERM, a higher star for- 
mation efficiency, f,/edc ~ 16 is required to fit the data at 
z = 6.56 for the LRM. As a result the SFR of LAEs in the 
LRM are increased by the same amount, ranging from 11 
to 197M0 yr~^. As in the ERM, we use the same value of 
fesc,a ~ 0.3 (increasing by 40% for larger halo masses) for 
z — 5.7 and 6.56, but the data at z = 4.5 again require fa 
to decrease by a factor of 4. 

A comparison between the Lya transmissivity, Ta, for 
the two reionization models considered is shown in Fig|4]for 
z — 6.56. In both cases the transmissivity increases towards 
more massive halos because of their generally larger SFR; 
also, at a given halo mass, Ta varies from 0.36 to 0.51 for 
the ERM, while it varies from 0.01 to 0.26 for the LRM i.e. 
it is considerably smaller for the LRM. In the LRM, small 
LAEs are characterized by a lower Ta with respect to larger 
ones relative to ERM. This is because even though the SFR 
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Figure 4. Lya transmissivity as a function of tlie LAE dark 
matter halo mass at ^ = 6.56 for the LRM (solid line) and ERM 
(dashed). 



are higher than in the ERM, the smaller LAEs are not able 
to build large enough HII regions; as a result, their Lya line 
is much more damped as compared to that for the larger 
LAEs. 

In conclusion, the LF data seem to require a strong 
increase of the SFR from z — 5.7 to 6.56 in the LRM to 
fit the observed LPs while a SFR that smoothly decreases 
with increasing redshift fits the observations for the ERM. 
Looking at the general trend, one finds that SFR densities 
decrease with increasing redshift. Hence, we find that the LF 
data favors the reionization scenario described by the ERM, 
i.e. a highly ionized (~ 3 x 10"*) Universe aX z — 6.56. The 
Best fit parameter values for the ERM are shown in Tab. 2. 

A caveat is that this analysis has been done for isolated 
emitters. As shown by McQuinn et al. 2007, clustering sig- 
nificantly increases the amount of Lya luminosity that can 
be transmitted by an emitter by adding a boost term to 
the background ionization rate. We find that such a boost 
factor of ~ 100 boosts the luminosity transmitted by the 
LAEs ai z — 6.56 with xhi = 0.15 significantly and in that 
case, the LRM can be fit by the same parameters {fesc,a, 
ft/^dc) as the ERM. However, an estimate of the boost in 
the background requires an accurate understanding of the 
radial dependence of the clustering and the contribution of 
each emitter to the boost. We then defer this analysis to fur- 
ther papers where we would use LSS simulations to fix these 
quantities. Hence, we can not rule out the LRM completely 
till clustering is included and better measurements of SFR 
densities at z > 6.56 are obtained. 



4.3 UV luminosity function 

Shimasaku et al. (2006) transformed the z' band magnitude 
from the photometric sample of 89 LAE candidates into the 
far UV continuum at the rest frame. The UV LF was calcu- 



lated by dividing the number of LAEs in each 0.5 magnitude 
bin by the effective volume corresponding to the FWHM of 
the bandpass filter used (NB816). Objects fainter than the 
2(7 limiting magnitude (27.04 mag) in the z' band were not 
included in calculating the UV LF and this corresponds to 
the vertical line at Muv = —19.58 in Fig[S] The authors 
mention that the apparent flattening at Muv > 20.5 might 
be due to the incompleteness in the measurement of the far 
UVLF. 

Kashikawa et al. (2006) used the same methodology 
mentioned above to derive the rest UV continuum from their 
photometric sample of 58 LAEs. Their LF measurements at 
magnitudes fainter than Muv = —20.24 (Sa) are uncertain 
due to the z' band magnitudes no longer being reliable be- 
yond this value. 

Both the above calculations have accounted for the de- 
tection completeness of the narrow band filters. They also 
find that cosmic variance is not severe for the UV LF. An 
important point to note is that the UV LFs a.t z = 5.7 and 
z = 6.5 are in very good agreement and show no evolution 
between these redshifts, which is in clear contrast to the Lya 
LF which shows a deficit of high luminosity LAEs at « = 6.5 
as mentioned before. 

We derive the specific continuum luminosities using 
STARBURST99, adopting Z = 0.05^q, an age of about 100 
Myr and a Kroupa IMF (details in Sec. 2.2). The contin- 
uum luminosity is then related to the SFR by 

Lc(1375A) = 2.13 X W'^°[M,/MQyr-'^]erg s"^l"\ 

However, using this conversion and the best fit parameter 
values of f*/edc for the ERM as mentioned in Sec l4.2l we 
find that the UV LFs for both redshifts lie above the ob- 
served ones. Hence, additional dust damping of the UV LF 
is required to match with the observation^J. We quantify 
this additional damping by introducing fc, the fraction of 
continuum photons that escape the LAE, unabsorbed by 
dust. Using a single value of fc for a specific redshift (see 
Tab[2}, across the entire mass range considered, we find a 
reasonably good agreement with the observed UV LF for 
the bright LAEs. However, the model fails to reproduce the 
bending of the UV LF observed for the low luminosity emit- 
ters. This could either be due to detection incompleteness in 
the observations or due to the lack of a physical effect such 
as a halo mass dependent escape fraction of UV photons. 
A simple prescription for the latter would be an increasing 
dust content with decreasing halo mass (due to a decrease 
in the ejection efficiency). However, other explanations such 
as SFRs that decrease with decreasing halo masses can not 
be ruled out with this model. A full exploration of possible 
effects will be carried out in further works using simulations. 
It is interesting to note that for this model, while at 
the highest redshift, continuum photons are less absorbed 
by dust as compared to the Lya photons, the trend reverses 
at lower redshifts. This could hint at dust whose inhomo- 
geneity /dumpiness evolves with redshift. However, robust 
estimates of the ages, metallicites, IMF and detailed stud- 
ies of dust distribution and its evolution inside LAEs are 
needed before such a strong claim can be made. 



We have taken the continuum luminosity value averaged over 
1250 to 1500 A, with the centre at 1375 A. 
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Table 2. Best fit parameter values for the ERM to fit both the 
Lycf LF and UV LF. For each redshift (col 1), we mention the 
halo mass range required (col 2), the SFR efficiency (col 3), the 
associated SFR (col 4), the efi'ective escape fraction of Lya pho- 
tons (col 5) and the escape fraction of continuum photons (col 
6). 



z 


MhWq] 


/*/.dc 


M.[M0yr-i] 


J esc, a 


/c 


4.5 


IQll.l- 


-12.5 


3.5 


6- 


-160 


- 0.075 


~ 0.045 


5.7 


IQlO.S- 


-12.3 


3.5 


3- 


- 103 


~0.3 


~0.25 


6.56 


1010.7- 


-12.0 


3.5 


2- 


-43 


~0.3 


~0.5 



-3.5 



-4.5 




-5.5 



-20 -20.5 -21 -21.5 
M„, (AB)fmag] 



23 



Figure 5. UV LAE Luminosity Function for the early rcioniza- 
tion model (ERM). Points represent the data at two different 
redshifts: z = 5.7 Shimasaku et al. (2006) (circles), z = 6.56 
Kashikawa et al. (2006) (triangles). Lines refer to model predic- 
tions at the same redshifts: z = 5.7 (dot-dashed), z = 6.56 (solid). 
The vertical dashed (dotted) lines represent the 2a (3cr) limiting 
magnitude for z = 5.7 {z = 6.56). 




Figure 6. Contribution of LAE to the cosmic SFR density evo- 
lut ion from our b est fit models. Points show the measurements 
bv lHopkinj ll2004h : the dashed (solid) line is the prediction from 
ERM (LRM). 



easy to interpret. As the dust formation timescale is about 
10 Myr, if the latter is copiously produced in supernova 
ejecta, as pointed out by several authors (Kozasa, Hasegawa 
& Nomoto 1991, Todini & Ferrara 2001, Schneider, Fer- 
rara & Salvaterra 2004, Bianchi & Schneider 2007) and re- 
cently confirmed by the extinction curves of high redshift 
quasars (Maiolino et al. 2004), supernova-produced dust 
would rapidly increase the opacity to both continuum and 
Lya photons, thus causing a rapid fading of the emitter. 

Second, the contribution of LAEs is about 8% of the 
cosmic star formation rate density at z = 5.7. Thus, either 
the duty cycle of the actively star forming phase in these 
objects is of the same order, or one has to admit that only 
a very small fraction (~ 1/12) of high redshift galaxies ex- 
perience this evolutionary phase. In the first case, the star 
formation duration would last about 8% of the Hubble time 
at 2 = 5.7, i.e. 72 Myr. 



4.4 Cosmic star formation rate density 

As a sanity check, using the parameters that best fit the 
data as discussed in Sec l4.2l we calculate the contribution 
of LAEs to the SFR densities at 2 = 4.5, 5.7 and 6.56. We 
compa re these with the SFR densities observed by iHopkinj 
1I2OOJ) (Table 2) for the common dust-correction case, the 
results for which are plotted in Fig|6] We find that for the 
best-fit parameters, the contribution of LAEs to the SFR 
density is redshift-dependent, being about 8% at z = 5.7 
with SFR in the range 3 < M*/M0yr"^ < 103, and even 
higher at z = 4.5, although the data present a large scatter 
at the latter epoch. 

Further, two points are worth noticing about the pre- 
dicted SFR density. First, the SFR density must increase 
strongly from z = 5.7 to 6.56 in the LRM case. Although 
not impossible, such behavior is certainly puzzling and not 



4.5 Lya equivalent width 

From our model it is easy to derive the intrinsic rest-frame 
Lya line equivalent widtqj. Since both the continuum and 
I/™' scale linearly with SFR in our model, the intrinsic EW 
distribution is a 5-function at EW^"^ « 13lA. 

From our model, the observed EW in the rest frame of 
the emitter is calculated as 



2.8 X 1O*2[M*/M0yr--Jerg s" 



-11 



^ We calculate the intrinsic rest-frame EW as 
LJ,"V[Lc(1375A)] where L™« 
and the specific continuum luminosity is given by Lc = 2.13 X 
lO*''[M,/M0yr~-'^]erg s~-'^A~-'^. Both the intrinsic Lya and con- 
tinuum luminosities have been derived using STARBURST99; we 
adopt Z = O.O5Z0 and a Kroupa IMF (for details, see Sec. 2.2). 
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Figure 7. Normalized distribution of the rest fram e EW f or LAEs 
Sut z = 4.5. Observed values from I Dawson et al.l (|2007|) (model 
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Figure 8. Normalized distribution of the rest frame EW for LAEs 
at ^ = 5.7. Observed values from Shimasaku et al. (2006) (model 
results) are shown by solid (dot-dashed) lines. 



EW = £]¥""'{!- f,sc)T^ { ^ 



(19) 



where /c quantifies the fraction of the continuum luminosity 
which escapes the emitter, unabsorbed by dust. The ratio 
fa/fc expresses the differential extinction of the Lya line 
with respect to continuum radiation due to dust grains. 

At z — 4.5, we have seen that we require a factor ~ 4 
suppression of the Lya line luminosity by dust, i.e. fesc,a ~ 



0.075. As dust affects also the continuum, and hence the 
EW, we need to estimate the value of fc (calculated at A = 
1375 A). We find that fc ~ 0.045 for the mean EW from 
our model (155 A) to be the same as the observed EW (155 
A). We then use the following relations to obtain the color 
extinction: 



AA(1375A) = -2.51og/c, 

Av 1Aa(1375A) 



E{B~V) 



Rv 



Rv 



(20) 
(21) 



where Rv ~ 3 and we have assumed a Galactic extinction 
curve. From these expressions we obtain E{B — V) = 0.28. 
The value of fc implies that the continuum is extincted 
about 1.6 times more heavily than the Lya line (assum- 
ing fc3c ~ so that /□, — 0.075). This is not inconceiv- 
able if LAE interstellar dust is inhomogenously distributed 
and/or clumped, as showed by Neufeld (1991). With these 
two values we then derive the predicted EW distribution and 
compare it with the Dawson et al. (2007) data in Fig[71 As 
mentioned before, for the best fit parameters to the LF at 
z = 4.5, fesc,a ~ 0.075 and Ta ~ 0.50. Note that fesc^a and 
particularly Ta depend on the LAE luminosity/mass and 
increase by about 45% and 20% respectively towards higher 
masses. 

The predicted EWs are concentrated in a range, 
114 A < EW < 201 A (mean=155 A), whereas the ob- 
served distribution is considerably wider, spanning the range 
6 — 650A with a mean of 155A. As explained above, the 
spread of the predicted EW distribution arises only from 
the corresponding spread of SFR (6-160 M0yr~^) required 
in order to match the LF at z — 4.5, via the dependence of 
Tc on the SFR. 

Calculating the rest frame EWs is easier a.t z — 5.7 
since we have an estimate of fc from the UV luminosity 
function as mentioned in TaljS] We calculate the EWs using 
fcscc ~ 0.3, fc ^ 0.25 and T^ « 0.37. As for z = 4.5, fcsc,c 
and Ta depend on the halo mass and increase by 40% and 
45% respectively towards higher masses. The calculations 
then yield EWs that range between 56-127 A. The mean 
from our model (~ 92.3 A) is much less than the mean 
value of 120 A, observed by Shimasaku et al., 2006. 

The narrow range (z=4.5) and lower mean (at z=5.7) of 
EWs calculated from our model can easily be explained by 
the fact that our model does not include inflows/outfiows, 
assumes an age of about 100 Myr for all the emitters and 
a metallicity which is 1/20 of the solar value. In reality, a 
larger spread would be expected from the addition of phys- 
ical effects lacking in this model, such as (i) gas kinematics 
(inflow/outflow); (ii) variations of the IMF, metallicity, and 
stellar populations (including PopIII stars), and (iii) young 
stellar ages. 

While infiows erase the red part of the Lya line, thereby 
reducing the EW, outflows shift the line centre redwards, 
helping more of it to escape. Outflows can also add a bump 
to the red part of the line due to backscattering of Lya 
photons, as shown by Verhamme, Schaerer & Maselli, 2006. 
A top heavy IMF produces more H i ionizing photons, as 
does decreasing the metallicity. Hence, both these effects 
increase the EW. Further, for very young emitters (~ 10 
Myr), the continuum is much less than the continuum at 
100 Myr and so, the EW would be much larger for younger 
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Figure 9. Weighted skewness of the observed Lyg line for dif- 
ferent models. The asterisks are the data from iKashikawa et alj 
1I2OO6I '). The dashed (solid) line correspond to the best fit ERM 
(LRM) at 2 = 6.56. 



emitters. All these effects need a much more dedicated study, 
which we defer to future work. 



4.6 Line profile asymmetries 

Additional constraints on the model can come from the in- 
formation embedded in the observed line profiles, as for 
example the line profile asymmetry. This can be suitably 
quantifie d by the weighted skew ness parameter, Sw, intro- 
duced bv lKashikawa et al.l l|2006r ). which we calculate for the 
best fit parameter values for the ERM mentioned above. We 
adopt the following definition for such a quantity: 



Sw = S/\X = SyXiQ^r — Ai 



(22) 



where Aio,r (-^10, b) is the wavelength redward (blueward) of 
the Lyo line where the flux falls to 10% of the peak value. 
In addition, we have that 



I = 



S = 



/ 

I 

Ja3 



(23) 

(24) 
(25) 
(26) 



where fi is the line flux in the wavelength pixel i whose 
coordinate is Xi, and the summations are performed over 
the pixels covered by the Lya line. On general grounds one 
would expect that the observed Lya line shape would be 
more symmetric (i.e. low Sw) in reionization models char- 
acterized by a lower value of xhi- However, given the above 
definition, just the opposite is true. In fact, for any reason- 



able value of the relevant parameters (see Fig.l) the blue 
part of the line is heavily absorbed, thus yielding a high 
value of Sw; as xhi is increased, also the long- wavelength 
part of the line is affected by the red damping wing, making 
the line more symmetric around the peak. 

The predicted trend of Sw with the observed Lya lumi- 
nosity ai z — 6.56 is reported in Fig[9] for the parameters of 
the ERM and LRM that best fit the LF data (discussed in 
the previous Section). For both models, the weighted skew- 
ness of the line increases for more luminous objects; however, 
such dependence is steeper for the ERM than for the LRM. 
In general, though, the two reionization scenarios predict 
Sw values in the range 10-17. The data from Kashikawa et 
al. (2006) spans the somewhat larger range 3-17, with many 
of the data points lying around Sw = 5. Given the paucity of 
the observed points and the large errors associated to them, 
it is probably premature to draw any strong claim from these 
results. However, given the constant increase in the amount 
and quality of LAE data, it is quite possible that the line 
skewness could represent a very interesting tool to constrain 
reionization models in the near future. It has to be noted 
that the data show a large scatter of Sw at a given value 
of La, perhaps indicating that local conditions, including 
gas infall/outffow, density inhomogeneities and interaction 
of the Lya line with the interstellar medium of the galaxy, 
might play a dominant role. These can only be investigated 
in a statistically meaningful manner via high-resolution nu- 
merical simulations to which we will defer a forthcoming 
study. 

From the theoretical point of view it is instructive to 
summarize the response of the skewness to different physical 
conditions. As we have seen from Fig[9l Sw increases with 
La (or, equivalently with SFR); this is true for any fixed 
value of Xhi- This is because as the SFR increases, more of 
the Lya line escapes forcing Sw to increase as a result of the 
larger value of AA. Further, the long -wavelength part of the 
observed Lya line begins to fiatten with increasing xhi due 
to attenuation by the red damping wing. Hence, AA varies 
slower with La (SFR) for high xhi (LRM) as compared to 
lower values (ERM); this makes the slope of Sw steeper for 
the ERM. 

A more general view of the dependence of Sw on La 
(hence on SFR) and xhi is shown in Fig. 1101 The plot has 
been obtained by dividing the observed Lya luminosity into 
bins and averaging the weighted skewness over the number 
of LAEs in each bin at a given value of xhi- The regions 
with weighted skewness values equal to zero represent a lack 
of LAEs in that bin. 

The most intriguing feature of Fig llUI is a clear anti- 
correlation between La and xhi- Given the range of SFR 
considered (M* = 2.7 - 197M0yr"^), LAEs populate pro- 
gressively fainter Lya luminosity bins as the IGM becomes 
more neutral. Notice that relatively luminous objects {La ~ 
lO*^'^ erg s"^) would not be detected ff xhi ~ 0.25. Within 
the range in which these objects are visible, the most lu- 
minous objects always show the largest Sw at fixed Xhi', 
however, such maximum value is also seen to increase with 
decreasing Xhi- 

The model does not include important effects such as in- 
ffows/outflows and interaction of the Lya photons with the 
ISM, which will definitely leave an imprint on the Sw and 
hence, weaken the La-Xni anti-correlation. However, a com- 
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Figure 10. Dependence of Sw (values are color-coded by the bar on the right) on xm and La at z = 6.56 for a set of LAEs with SFR 
in the range predicted by the two reionization models, i.e. M* = 2.7 — 197MQyr~'^. 



posite Lya line, built from the observations of a sufficiently 
large number of LAEs might show such an anti-correlation. 
In future work, we will endeavor to include these effects, 
hence obtaining an estimate of how large a sample might be 
sufficient for this purpose. 



5 DISCUSSION 

Starting from a simple yet physical model of galaxy for- 
mation within dark matter halos coupled with a popula- 
tion synthesis code, we have derived the intrinsic Lya lu- 
minosity for a LAE. We then compute the volume of the 
ionized region built by the source and the density profile 
of the neutral hydrogen within it to obtain the damping 
of the emitted Lya line caused by the Gaussian core and 
Lorentzian wings. Using this semi-analytic model, we have 
first explored the physical dependence of the observed Lya 
line profile on various free parameters such as the LAE star 
formation rate, ionized region radius, the effective escape 
fraction of Lya photons and the global IGM neutral hydro- 
gen fraction. Among other things, we pointed out the in- 



teresting fact that the observed Lya luminosity, La, peaks 
at a value of fesc ~ 0.5 if the gas is substantially neutral 
(Xhi ~ 0.15); for that value, the contribution of fesc to the 
intrinsic Lya luminosity and the Stromgren sphere balance 
each other. 

By considering two physically motivated scenarios in 
which reionization occurs either early (ERM, Zi ~ 7) or 
late (LRM, Zi ~ 6) we have fixed the global value of the 
IGM neutral fraction, xhi (thus leaving the star formation 
efficiency and the effective escape fraction of Lya photons as 
the only free parameters), and obtained both the observed 
Lya line profile and the Lya luminosity function. Finally 
we have compared these predictions with available data at 
various redshifts. 

Using this procedure we have been able to fit the LPs 
obs erved by Dawson et al. (2007), Shimasaku et al. (2006) 
and lKashikawa et al.l (|2006r ) at z = 4.5, 5.7 and 6.56 respec- 
tively for the ERM. According to this model, no redshift evo- 
lution or mass dependence of the star formation efficiency 
is required. On the contrary, the LRM requires an increase 
of a factor 16/3.5 ~ 4.5 in the SFR efficiency from z = 5.7 
to 6.56. Although not inconceivable, such an upturn of the 
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star formation efficiency seems puzzling and at odds with the 
observed cosmic star formation rate density. On this basis, 
we are more inclined to support the ERM. In addition, we 
find that the evolution of the observed luminosity function 
from z — 5.7 to 6.5 does not imply that we are scratching 
the reionization surface yet. Rather, the LF evolution can 
be explained solely by the evolution of the underlying dark 
matter halo mass function between these redshifts, as has 
previously been discussed by Dijkstra, Wyithe & Haiman, 
2007. One would however, require more observations of the 
SFR density and information regarding the boost added to 
the ionizing background due to clustering at 2 ~ 6.5 to com- 
pletely rule out the LRM. 

A reasonable fit to the data at z — 5.7 and z = 6.56 is 
obtained for a single value of fesc,a ~ 0.3 (although a good 
fit is obtained by allowing for a (40%) increase of fesc,a to- 
wards larger masses). The data at 2 = 4.5 instead pose a 
challenge to the model, as outlined earlier, because the ob- 
served number density of luminous objects is lower than that 
predicted by the evolution of the theoretical LF at higher 
redshifts. This could imply an increase in the overall dust 
content of LAEs at this redshift which would lead to absorp- 
tion of Lya photons within the emitter. Such a clumpy dust 
component is also suggested by by the large EWs observed 
at 2 = 4.5. 

We obtain the UV LF at 2 = 5.7 and 6.5 for the best 
fit values of f*/(.dc for the ERM. We find that additional 
damping of the UV luminosity is needed to match the pre- 
dictions with the observations and quantify this by /c, the 
escape fraction of continuum photons. A single value of /c 
for a given redshift is enough to match the high luminos- 
ity end of the UV LF but does not produce the bending 
required at the faint end. However, this bending might just 
be the result of detection incompleteness of the sample. An 
interesting result is that for the given IMF and metallicity, 
while at higher redshifts (2 = 6.5), the continuum photons 
are less absorbed by dust as compared to Lya photons, this 
trend reverses at lower redshifts(2 = 4.5, 5.7). This could 
be explained by imhomogenously distributed/clumped dust. 
However, the IMF, ages and metallicies of the emitters must 
be fixed robustly using simulations and infall must be in- 
cluded in the model before such a strong claim can be made. 

Using the best fit value of / * /edc and escape fraction 
of Lya photons (continuum photons) obtained for the ERM 
from the Lya (UV) LF, we calculate the expected EWs at 
2 = 5.7 and find that the mean (~ 92 A) is much less than 
the observed value of 120 A. At 2 = 4.5, since there are 
no observations of the UV LF at present, we calculate the 
escape fraction of continuum photons required to match the 
EW mean from the model to the data. A dust extinction 
oi E(B — V) ?5i 0.28 brings the predicted mean Lya EWs 
(« 155 A) in very good agreement to the observed mean 
(« 155 A). This value of dust extinction is reasonable when 
compared to the observational upper limit of E{B — V) ~ 0.4 
(Lai et al., 2007). However, additional efi'ects which vary on 
a galaxy to galaxy basis, such as outflows/inflows or peculiar 
stellar populations are required to account for the spread of 
EW seen in the data. 

The contribution of LAEs to the cosmic SFR density 
is small, amounting to roughly 8% sd, z — 5.7. Thus either 
the duty cycle of the actively star forming phase in these 
objects is of the same order, or one has to admit that only 



about one-twelfth of high redshift galaxies experience this 
evolutionary phase. 

Additional useful information can be extracted from the 
line profile by using indicators like the line weighted skew- 
ness and equivalent width. The results presented here (Sec- 
tions 4.4-4.5) must be considered as very preliminary for 
several reasons. First the available data on both Sw and 
EW are very scarce and of relatively poor statistical qual- 
ity, as they are very difficult to obtain even from the best 
current observations. Second, our model contains a number 
of simplifications which make the comparison only meaning- 
ful at a basic level. Nevertheless, it is encouraging that the 
model results are broadly in agreement with the data, at 
least for what concerns mean values. 

The models presented in this study do not include feed- 
back processes related to the energy injection by supernovae. 
As pointed out by Santos (2004) and Ificv ot al. (2007), pe- 
culiar gas motions might affect the line profile consider- 
ably: while infiows of gas erase the Lya line, galaxy scale 
outflows produced by supernova (or AGN) feedback enable 
more of the Lya to escape. Yet, our models are able to 
fit the LF evolution in the redshift range 4.5 < 2 < 6.56 
quite well. This might indicate that the effect of feedback 
might be similar for all the emitters in the halo mass range 
(Afh = 10"''^"^^"AfQ at 2 = 6.56). Obviously, a more firm 
statement can be made only after a proper inclusion of pe- 
culiar motions. The backside of this is that, if infiows are 
taken into account, o ne cannot cons train the value of xhi 
robustly, as noted by ISantosI (|200J) . These effects can be 
properly taken care of by using high-resolution numerical 
simulations, which we plan to use in future works. 

In the same spirit, a full study of the problem should 
also include the effects of IGM density and temperature in- 
homogeneities, precise values of the local metallicity and 
star formation rates, and require information about the spa- 
tial clustering of the emitters. These values will be fixed in 
subsequent papers using the results of the simulation by 
iTornatore. Ferrara fc Schneide"3 (|2007l '). 

An even more cumbersome ingredient is represented by 
dust. Dust grains act as sinks of Lya and continuum pho- 
tons, thus depressing the Lya line luminosity but possibly 
boosting the line EW as explained throughout the paper. 
To what extent and on what timescales LAEs become dust- 
polluted (and possibly enshrouded) remains a question to 
which both theory and observations can provide only coarse 
answers at this time. Our conclusions hint at the need for 
dust in order to explain the evolution of the LF toward the 
lowest redshift. However, lacking a precise knowledge of the 
mass dependence of the dust-to-gas ratios in high redshift 
galaxies and a deep understanding of the dust formation 
processes/sources, developing a fully consistent theory will 
keep us busy for many years to come. 
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